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Abstract 

We discuss the description of a many-body nuclear system using Hamilto- 
nians that contain the nucleon relativistic kinetic energy and potentials with 
relativistic corrections. Through the Foldy-Wouthuysen transformation, the 
field theoretical problem of interacting nucleons and mesons is mapped to an 
equivalent one in terms of relativistic potentials, which are then expanded at 
some order in 1/mjv. The formalism is applied to the Hartree problem in 
nuclear matter, showing how the results of the relativistic mean field theory 
can be recovered over a wide range of densities. 
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I. INTRODUCTION 



In recent years considerable efforts have been devoted to the development of relativistic 
nuclear models, whose basis is usually a relativistic Lagrangian containing nucleons and 
mesons . On the other hand, the traditional concept of interparticle potential has proved 
to be quite useful in nuclear physics and rather sophisticated techniques of calculation have 
been developed within the hamiltonian formalism. 

It is clear that it would be useful to be able to reformulate the relativistic problem of 
interacting Dirac particles in terms of Pauli spinors and potentials, since this would allow 
one to employ the non-relativistic computational techniques. Furthermore, the high energies 
and high momentum transfers accessible at facilities such as Cebaf require calculations to 
be performed in regimes where the relativistic kinematics is certainly important: then, one 
may ask the question whether genuine relativistic dynamical effects are present and to which 
extent they may reconducted to a potential treatment. 

A possible way of connecting the two approaches relies on Hamiltonians containing the 
nucleon relativistic kinetic energy and a potential obtained by expanding at some order in 
l/rriN the nucleon-nucleon (NN) scattering amplitude. For instance, in Ref. [|[ the quasielas- 
tic charge response of electron scattering has been evaluated in the relativistic Fermi gas 
model using the Bonn potential expanded up to A similar approach has been followed 

in Ref. ||^, using a relativistic Hamiltonian in variational Monte Carlo calculations. 

The aim of the present paper is to show how a systematic expansion of the nucleon- 
nucleon interaction in a many-body nuclear system can be constructed, using the well- 
known Foldy-Wouthuysen transformation. The many-body problem, expressed in terms 
of relativistic Green's functions and meson-nucleon vertices, is transformed to an exactly 
equivalent one, in terms of modified Green's functions and vertices. After projection on the 
positive energy states and expansion in powers of I/ttitv one recovers the familiar relativistic 
potentials. 

This procedure gives a systematic way of dealing with l/rriN expansions in a many-body 
context: it allows one to address, for instance, the issue of self consistency, i. e. the fact that 
the Dirac spinors depend on the nucleon mass, which is modified by the medium. In terms of 
relativistic Hamiltonians, this means that the potential, which affects the effective nucleon 
mass via the self-energy, should in turn depend on it. In a fully relativistic context this 
problem is rather simple to solve, whereas in a non-relativistic framework it is ill-defined. 

In the following we shall not be concerned with antinucleons degrees of freedom: although 
this is an interesting issue in itself, we would like to maintain our initial commitment of trying 
to understand the nuclear phenomena including "simple" relativistic effects, such as kinetic 
energies and momentum-dependent potentials. Any failure in this program might of course 
open the way to more interesting relativistic effects, like, e. g., vacuum polarization. 

The reader should also be aware of a potential source of inconsistency connected to 
retardation effects. In fact, the parameters of empirical NN potentials, such as the Bonn one, 
have usually been fixed assuming a static interaction. Although the inclusion of retardation 
in the formalism developed in the present paper is rather straightforward, when dealing 
with an empirical NN potential this may not be the most appropriate course to follow. 
Since in the following we have applied the formalism only to Hartree calculations, which are 
not affected by this problem, we have considered, for simplicity, static propagators for the 
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meson fields. The issue will have to be reconsidered when dealing, e. g., with random phase 
approximation calculations. 

The paper is organized as follows: in Sect. 0, the two most commonly used schemes 
for defining a relativistic potential are briefly sketched, employing a and uj exchange as 
examples. In Sect. |T^, the Foldy-Wouthuysen expansion is introduced in the nuclear many- 
body system, using nuclear matter for illustration, whereas in Sect. |^the Hartree problem 
is solved, comparing the results obtained with the relativistic Hamiltonians to the exact 
ones. Finally, in the last Section we present our conclusions. 



II. RELATIVISTIC NN POTENTIALS 

In the literature one can find essentially two ways to derive a two-body potential ex- 
panded up to a given order in 1/mjv. 

In the first case one expands the relativistic scattering amplitude and interprets it as the 
matrix element of the Fourier transform of the potential (Born term) between Pauli reduced 
spin wave functions. We shall refer to it as "Breit reduction" 

An alternative procedure starts from the three-dimensional Blankenbecler-Sugar (BbS) 
reduction of the Bethe-Salpeter (BS) equation, which yields an amplitude satisfying elastic 
unitarity 

Both derivations are briefly sketched in the two following subsections. Here, to set the 
notation, let us define the momentum transfer q, the relative momentum in the initial state p 
and the total momentum P for the scattering from a two nucleon initial state with momenta 
fci and k2 to a final state with momenta k'^ and k'2. 

q = ki-k[ = k2- fca (2.1a) 
p=^{k,-k2) (2.1b) 
P = ki + k2 = k[ + k'2 . (2.1c) 

We will be considering the amplitudes for the exchange of a scalar particle of mass 
(described by the scalar field 0) and of a vector particle of mass rriy (described by the vector 
field v^), whose couplings to the nucleon are given by 

= Qs^^^ (2.2a) 
= g.'^Y'^v^ , (2.2b) 

where \1/ is the nucleon field and 'j^ are the Dirac matrices. 

A. Breit reduction 

The amplitudes for the exchange of a scalar meson, Mg, and of a vector meson, are 
given by: 
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where u(h) is the Dirac spinor, 



u{k!^-i^u{kx)u{k!^-i^u{k2) 



(2.3a) 
(2.3b) 



u{k) 



2m N 



X 



(2.4) 



whereas x is the Pauh spinor and Ek = \J k? + m?^. 

When expanded up to order k^/mj^f, the normahzed Dirac spinor turns out to be 



u{k) 



5m 
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crk 

2mjv 



X 



+ 



mN, 



(2.5) 



Inserting this expression in (|2.3|) and using the standard relation between the Born 
amphtude and the spin matrix elements of the potential in momentum space, M = 
x'lx'lViq,p, P)X2Xu one gets 
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[2.6a) 



(2.6b) 



where the superscript BR stands for "Breit reduced" and Su is the standard tensor operator 
5'i2 = 3(cri ■ q){(T2 ■ q)/q^ - cr^ ■ (72. 



B. Minimal relativity 

A natural framework for relativistic two-nucleon potential scattering is given by the BS 
equation. This is a four-dimensional integral equation for the scattering amplitude, which is 
rather hard to solve even in the ladder approximation. Hence, many approximation schemes 
have been derived in order to obtain from it a more tractable three-dimensional, covariant 
equation, based on the requirement of relativistic elastic unitarity 

A popular reduction procedure leads to the BbS equation, which, in the centre-of-mass 
(cm.) frame, reads 
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M{kf,ki) = V{kf,ki) 



mi, 



(2), 



(2.7) 



where V is the bare two-body potential, A^'' are the positive energy projection operators 
and fcj, k and kf are the initial, intermediate and final cm. momenta, respectively. 
If one defines 



M{kf,ki 
V{kf,k, 



irriN 
Eh, 



M{kf,k, 



V{kf,k, 



hi 



E, 



ki 



(2.8a) 
(2.8b) 



Eq. 



becomes identical to the non-relativistic Lippmann-Schwinger equation 

A^^\k)A^^\-k) 



M{kf, ki) = V{kf, ki) + TUN / dkV{kf, k) 



kf — + ie 



-M{k,ki 



(2.9) 



Prescription ( |2.8| ) is usually known as "minimal relativity" It can be generalized to a 
generic frame by introducing a factor (mjvZ-E)* for each spinor in the amplitude, i.e. 



mjv mN TTlN TTlN 

Ek-i Ek'^ Ek^ Ek'^ _ 



^ 1 

4 



V{q,p,P) 



(2.10) 



where MR stands for "minimal relativity" . 

Expanding again up to the order k'^/m'lf one gets, for the exchange of scalar and vector 
mesons, the following expressions 



yf^(g,p,p) 



V^\q,p,P) 



9l 



+ ml 



{2p-qf 



5m 



N 



5m 



N 



5m 



N 



—ip ■ 



[q X (cTi + cTaJ 

Am% 



iP 



[q X (cTi - 0-2)] 



5m 



N 



1 - 



9I 

+ ml 



■im 



N 



im 



N 



5m 



N 



Am% 

Cl ■ 0-2 + ^012 



6m^ 



12m2^ 



(2.11a) 



:2.11b) 



Eq. (|2.9| ) has been used in deriving the relativistic Bonn potential 0. Hence, the non- 
relativistic expansion of this potential, useful in nuclear structure calculations, is given, for 
a and uj exchange, by Eq. ( p.ll| ). 
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III. FOLDY-WOUTHUYSEN EXPANSION 



In this section we would like to rederive the previous expansions using the language of 
the Green's functions and within the scheme devised by Foldy and Wouthuysen ||T0| (FW) 
to decouple the large and small components in the relativistic wave function. 



Let us start by introducing the FW unitary transformation |jTT]JT2 

/ Ek + rriN ( ^ , 7 ■ fc \ 



T{k) 



2Ek 



U 



Ek + rriN, 



(3.1) 



and the associated FW nucleonic field operator, 



FW 



(3.2) 



(having used the property ■y^T'^'-y'^ = T). 

Using (|3.2|) one can define a FW Green's function 

iG^-^ix^y) ^ (vi/olT [vi/^W(x)#^(i/)] IvPo) 
= TiG{x, y)T . 

In momentum space, and for the case of the relativistic Fermi gas, the latter reads 
Gl'^ik, ko) = T{k)Go{k,ko)T{k) 



(3.3) 



d{k - kp) 



^{kp - k) 



Ek + ie ko- Ek 



IE 



1 



■ fco + Ek - is 



(3.4) 



where P± = (1 ± 7°)/2 are operators that project out the large/small components in the 
wave function and kp is the Fermi momentum. 

It is evident that one gets two pieces, acting separately on the large and small compo- 
nents, since the transformation (|3.1| ) has washed out all the operators inducing a mixing 
(i.e. the 7i's). 

It is also clear from the definition (|3.3|) that one can always redefine any given Feynman 
diagram in terms of and a transformed interaction (see Fig. |l|). 
For a typical meson-exchange interaction one has, for instance. 



V = g^T^D{q)T2 



(3.5) 



where g is the coupling constant, D{q) the meson propagator and Fj the vertex operators 
(e.g., U or 7^^). 

Then, the FW interaction would read 



^Fw ^ g^r\k',)TrT\k^)D{q)T\k',)T2T\k, 



(3.6) 



We make now a further simplification, neglecting in the following the contribution of antin- 
ucleons: 
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FIG. 1. Scattering amplitude in terms of bare and FW propagators and potentials. 



ko — Ek + is ko — Ek — is 



(3.7) 



It is a simple matter to move the operators P+ from the Green's functions to the potential, 
defining 



y.^Y = 9'P+T\k',)V,T\k,)P^D{q)P^T\k',)T,T\k,)P^ 



(3.8) 



The many-body problem is then formally equivalent to the standard non-relativistic one, 
in terms of a potential V^^^ and the Green's function 



d{k - kp) 



d{kp - k) 



ko — Ek + is ko — Ek — is 



(3.9) 



which is identical to the non-relativistic one, apart from the energy-momentum relation. 

For the case of scalar and vector meson exchange, the FW reduced potentials have the 
following form: 



T/FW 
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. (3.10b) 



It is then a straightforward matter to check, by expanding these expressions up to second 
order in the inverse nucleon mass, that one gets the potentials of Eq. ( |2.6| ). On the 
other hand, if one performs the expansions dropping in (|3.10| ) the big square root containing 
kinematical factors {Ek + mN)/'2'Ek, then one recovers the potentials of Eq. 

The kinematical factors of the minimal relativity prescription ( p.8|) have been introduced 
in order to reduce the BS equation to a non-relativistic Lippmann-Schwinger equation: 
hence, it would appear natural to use the non-relativistic expression for the kinetic energy 
in doing calculations with the potentials. Note, however, that in the literature both 
the non-relativistic [[7|,p!3| and the relativistic [0,0 kinetic energies have been employed. 



IV. RELATIVISTIC HARTREE CALCULATIONS 

We have seen that the FW transformation leads to a new Green's function, in which 
the small and the large components are decoupled. Simultaneously, we define the potentials 
according to the prescription (|3.8D. In this way we are able to operate within a framework 



that clearly resembles the non-relativistic one: the FW Green's function, projected upon 
the large components, has indeed the same structure of the familiar non-relativistic Green's 
function (see ( p.9| )), apart from the energy of the fermion, which is now fully relativistic. On 
the other hand, the potentials obtained from a non-relativistic expansion of V^^^ reproduce 
the ones obtained following other procedures (for example, see |p|,Pjr3|). 

Let us now go further and apply the FW transformation to the relativistic Hartree 
approximation for the nucleon Green's function, assuming again that only scalar and vector 
isoscalar mesons are exchanged among the nucleons. This is the basis of many relativistic 
nuclear structure calculations. We will compare the results in the present approach with 
those obtained in the mean field theory |^ (MET), i.e. the Hartree approximation in which 
antinucleon contributions to the self-energies are neglected. 

In the MET the vector field shifts the energy of the nucleon, while the scalar field dresses 
the mass. It is worth stressing that this property follows directly from the Lorentz structure 
of the self-energy: the scalar self-energy, which is a scalar under Lorentz transformations, 
transforms as the mass, while the vector self-energy transforms as the time component of 
a four-vector. An important related issue in a Hartree calculation is connected to self- 
consistency. Indeed, in a relativistic framework the spinors depend on the nucleon mass, 
which is modified in the medium and should then be calculated self-consistently. This is 
rather straightforward in the MET, because of the above mentioned Lorentz structure of 
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the self-energy. In a non-relativistic calculation, based on an effective potential expanded at 
some order in A;^/m^, the dependence on the nucleon mass has been shifted to the potential: 
however, in this case one can no longer rely on general symmetry arguments and the FW 
expansion provides a systematic way to deal with self-consistency. 

Let us start by writing the Dyson equation for the inverse Hartree Green's function 

Gn{k, ko)-' = Go{k, ko)-' - Sh , (4.1) 

where Go{k, k^)'^ = 7°A;o — 7 • fc — mjy is the inverse of the free Green's function, whereas 
Sh is the Hartree self-energy due to the exchange of a and u, i. e. 

Sh = — 7oSh 

= ^4 / yf^Tr [Gn{k, ko)] e*^ - f yf^ToTr [^'Gn{k, ko)] e^'^^ . (4.2) 

The solution of (|4.1| ) is clearly: 

Gnik, ko)-' = 7% - 7 . fc - m^, (4.3) 



with 

ko = ko + (4.4a) 
mlf = mN + (4.4b) 

being the nucleon energy shifted by the vector field and the nucleon mass dressed by the 
scalar field, respectively. These relations, together with Eq. ([4.2|) , define the self-consistent 
Hartree solution in the MFT. 

Applying the FW procedure, one can introduce in the same way a Dyson equation for 
the FW Hartree Green's function: 

Gl^ik, ko)-' = T^{k)T{k)Gl'^{k, ko)-'T{k)Ti{k) - T^{k)^nT^{k) , (4.5) 

where Gq^(A;, ko)-' = 7°/co ~ -E'fc is the inverse of the free FW Green's function, while T}i{k) 
is the operator (|3.1|), evaluated for a nucleon with dressed mass m^. The solution to this 
equation can again be written as 

Gl'^ik,ko)-'=^%-E:, (4.6) 



with = yk'^ + and ko and still given by Eq. ( [4.4|) ; that is, at this stage one 
recovers the results of MFT. 

Now, in order to make contact with the effective potentials of the previous Sections, one 
has to project away the small components and to express the self-energies in terms of the 
potentials. Starting again from Eq. (|4.5| ), the left hand side becomes 



P+Gl'^ik, ko)-'P+ = P+{ko - El) ^ P+Gn(+){k, ko)-' , (4.7) 
whereas the last term on the right hand side can be rewritten as 
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In the previous equation we have introduced for brevity V^J^-^{k,p) = V^{q 
p)/2, k + p). Then, Eq. (|4.5| ) becomes 



ko - El = {ko + S^) - 



k^ + m*j^{mN + Efj) 
El 



(4.8) 
0,(fc- 

(4.9) 



where we have gathered separately the terms shifting the energy and the terms dressing the 
mass. The self-consistency condition is still given by Eq. ( 4.4b| ), as one can easily check. 
However, we can now reformulate it in terms of the FW potentials, namely 



rrij^ = rriN — i- 



El 



m 



N 



(2^ 



Tr 



(4.10) 



If we expand the potential up to 0(l/m^), then we obtain an approximate self-consistency 
equation for the effective mass: using the a potential of Eq. ( |3.10a| ), that is, after expansion, 
the V;^^ potential of Eq. (giGa]) , we get 



m 



N 



niN - 1 



3k% 



mi 



10m^2 



(4.11) 



where p = 2/c|,/37r^ is the density of nuclear matter. 

The exact expression ( |4.4b| ) and the approximate one ( [4.11D may be compared with the 



standard non-relativistic formula for m^: 
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rriN 
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1 + 



k dk 



(4.12) 



Note that from this formula one gets a constant effective mass in the Hartree approximation 
only when the interaction is quadratic in the non-local part, whereas ([4.11|) yields a constant 



m 



N 



at any order. Using the and the potentials of Sect. || one gets 
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1 + 



It is apparent that while the V^^^ potential complies with the relativistic symmetry require- 
ments (i. e., only a exchange dresses the mass), in the case of both a and u exchange 
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The effective mass as a function of the Fermi momentum for the MFT of Eq. ( [4.4b ) 



(solid), the approximation (4.11) (dashed) and the non-relativistic models (4.13a) and (4.13b) 



(dotted and dot-dashed, respectively). The upper panels correspond to the calculation neglecting 
self-consistency (see text), whereas in the lower ones the relativistic m^'s are self-consistent. The 
left hand panels correspond to the MFT choice of parameters and the right hand ones to the Bonn 
potential parameters. Note that in (a) and (b) the solid and dashed lines practically coincide. 



are contributing. We display in Fig. □ function oi kp for two choices of the couphng 

constants: one set is taken from the MFT (panels (a) and (c)), where they are fixed 
in order to reproduce the saturation point of nuclear matter, and the other from the Bonn 
potential (panels (b) and (d)), which is fitted to the NN scattering dataQ. The values of 
the parameters can be found in Table |. The largest difference is, of course, in the vector 
meson coupling, which is artificially reduced in the MFT in order to simulate the effect of 
short-range correlations. 

In the case of the two upper panels ((a) and (b)), Eqs. ( |4.4b| ) and ( |4.11| ) have been solved 
dropping the requirement of self-consistency, i. e. utilizing the bare mass rriN instead of 
in the right hand side of these equations. In the lower panels ((c) and (d)) the fully self- 
consistent solutions have been given. In the non-relativistic cases we have always applied 



Note that the Bonn potential contains also form factors, depending on the momentum transfer, 
which could be easily accomodated in the formalism. However, since in the Hartree self-energy 
{q = 0) they produce only a slight rescaling of the coupling constants, we shall neglect them in the 
following discussion, which has only illustrative purposes. 
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TABLE I. Masses and coupling constants of the a and uj mesons in the MFT [Q| and in the 
Bonn potential (BP) g. 





rris (MeV) 




(MeV) 


5^/4^ 


MFT 


550 


7.29 


783 


10.84 


BP 


550 


7.7823 


783 


20 



Eq. ( CT) . 

First of all, one should notice the remarkably good agreement between the exact expres- 
sion ( [4.4b|) for and the approximate one (|4.11|) (in particular, in the non-self-consistent 
instance the two curves are practically indistinguishable). The non-relativistic approxima- 
tion fails badly at the standard nuclear matter density and above. Note, however, that when 
the Fermi gas is applied to relatively light nuclei an effective lower value for hp should be em- 
ployed. For instance, in Ref. 0] an analysis of the charge response has been performed, using 
the expanded Bonn potential, for ^^C at kp = 225 MeV/c. In that case, the non-relativistic 
m*^ differs from the exact one by less than 3% (Fig. |^d). 

The dependence on the parameters is rather mild: while this is expected for the relativis- 
tic calculations and for V^^^, which depend only on the rather stable a coupling constant 
(see Table D, it turns out that also the dependence on g"^ in V^^^ is not dramatic. 

It is worth noticing that in spite of the fact that the effective mass stemming from V^^^ 
violates the relativistic symmetry requirements, the numerical values turn out to be rather 
close to the ones obtained from V^^^: the contribution of the a meson is reduced to make 
room for the u meson. We compare in Fig. ^ the effective masses (|4.13|) , showing also the 
individual contribution of a and u in the case of V^^. 

Coming back to the approximate relativistic expression ( [l.llj ), it is clear from Figs. ^ and 
^ that the expansion breaks above kp ^ 300 MeV/c. While this is enough for most nuclear 
matter calculations, it is worth pointing out that the range of validity can be increased by 
considering higher order contributions. Note that for momentum space calculations, such 
as those of Ref. 0] , the use of a potential expanded up to an arbitrary order does not pose 
any major problem. 



V. CONCLUSIONS 

In this paper we have presented a systematic way of dealing with l/rriN expansions in a 
relativistic many-body nuclear system. The formalism has been applied to the Hartree prob- 
lem, showing the importance of accounting, at the same order in the expansion, also for the 
condition of self-consistency. Treating the expanded potentials strictly non-relativistically 
can lead to substantial discrepancies with the relativistic calculation, depending on the den- 
sity. On the other hand, a correct treatment of the self-consistency condition considerably 
enlarges the range of densities where the relativistic potentials are applicable, already at 
order l/m%. 

A word of caution should be spent about the use of an empirical NN interaction, such as 
the Bonn potential: as we have seen, the minimal relativity prescription, which is used in 
deriving the potential, does not appear to be consistent with a relativistic dispersion relation 
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FIG. 3. The effective mass as a function of the Fermi momentum for the non-relativistic models 
( 4.13aD (solid) and ( 4.13b| ) (dot-dashed) using the Bonn potential set of parameters. The effective 
mass stemming from when only the a (dotted) or the lo (dashed) mesons are employed is also 
shown. 



for the kinetic energy and, furthermore, it generates contributions to the nucleon effective 
mass stemming from vector meson exchange. On the other hand, using the Bonn potential 
parameters with a "Breit reduced" potential would not be, strictly speaking, consistent. 
While this may not be a major problem in practice (at least in the case of m*^ the dependence 
on the choice of parameters is rather mild), a better option might turn out to be the one of 
fixing the parameters on some nuclear properties. 

The results presented here indicate that at the mean field level it is possible to map, 
with high accuracy, the many-body relativistic field theoretical problem to an equivalent 
one, expressed within the standard hamiltonian formulation. The next step would be, of 
course, to verify whether such a correspondence is valid also for correlations beyond the mean 
field, like, e. g., those described by the random phase approximation (RPA). Calculations 



based on both non-relativistic |]I6[ and relativistic |^ Hamiltonians have in fact shown the 
importance of the exchange terms in the RPA series, whereas fully relativistic calculations 
have always been performed keeping only the direct (ring) contributions. 
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